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O " 

O ■ Abstract 

^ ■ Recently constructed gravity solutions with Schrodinger symmetry provide a new example of 

AdS/CFT-type dualities for the type of non-relativistic field theories relevant to certain cold atom 
systems. In this paper we use the gravity side to calculate n-point correlation functions of scalar 
fields by reducing the computation to that in ordinary AdS space via a particular Fourier transform. 
We evaluate the relevant integrals for 3- and 4-point functions and show that the results are 
^ , consistent with the requirements of Schrodinger invariance, the implications of which we also work 
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out for general n-point functions. 
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I. INTRODUCTION 

Gauge/string duality is a central theme of modern research in string theory . As a 

strong/weak coupling duality it provides us with a wealth of data useful for studying inter- 
esting field theories at strong coupling where other methods fail. Its applications currently 
range from nuclear physics to plasma physics to tabletop condensed matter systems 4|. 

Until recently the study of gauge/gravity duality has been limited to relativistic field 
theories. However many non-relativistic field theories clearly play important roles in various 
physical systems. An example of a system with non-relativistic conformal symmetry is a 
system of fermions at unitarity which can be realized experimentally in certain cold atom 
systems [s!]. This theory is invariant under a non-relativistic conformal group called the 
Schrodinger group, and it is natural to wonder if there exists a gravity dual which would 
allow us to use the AdS/CFT correspondence to learn more about these systems. 

Recently a new exciting example of such duality has emerged. The simplest geome- 
tries with Schrodinger group isometry have been constructed in , leading to what is 
sometimes referred to as the AdS/cold atom correspondence or non-relativistic AdS/CFT. 
Excitingly, these geometries have been embedded in string theory in j^, y, Q . These papers 
have opened up a new avenue for studying non-relativistic conformal field theories (NR- 
CFTs) at strong coupling. See 
AdS/CFT. 

The holographic dictionary has been developed in [6|, 17| and the two-point function of 
a scalar field was recovered from a gravity calculation in g, 0]- Unlike more familiar ap- 
plications of AdS/CFT where the boundary is one dimension less than the bulk, here the 
boundary theory has two fewer dimensions than the bulk 



ll| for recent studies of various aspects of non-relativistic 



In this paper we consider hig. 
employ the trick introduced in 



ler point correlation functions in the gravity theory. We 



12| where the free massive Schrodinger equation was re- 



duced to a non-massive Klein-Gordon equation via Fourier transform with respect to the 
mass. Applying the same trick to the solution of the scalar wave equation in the gravity 
background for non-relativistic field theory, we reduce the wave equation to that in ordinary 
AdS space. The computation of correlation functions in the bulk then reduces to perform- 
ing a particular Fourier transform of AdS correlators. We explicitly evaluate the relevant 
integrals for three- and four-point functions and show that the results are indeed consistent 
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with the requirements of Schrodinger invariance, the imphcations of which we work out 



for general ra-point functions. The same integral for three-point function appeared in |12| . 
where these nonrelativistic three-point functions where identified with response functions in 
Martin-Siggia-Rose theory, but the integral has not been explicitly evaluated there. 

The paper is organized as follows. In section 2 we briefly review the necessary results 
from non-relativistic conformal field theories and discuss the constraints that Schrodinger 
symmetry places on n-point functions of primary operators. In section 3 we compute the 
bulk-to-bulk and bulk-to-boundary propagators for a scalar field. In section 4 we compute 
three- and four-point functions on the gravity side and compare them with the corresponding 
correlators in NRCFT. In appendix we explicitly evaluate the gravity bulk integrals which 
appear. 



Note Added. While this paper was in preparation the paper [13| appeared which also 
considers three-point function in non-relativistic AdS/CFT. 

II. CORRELATION FUNCTIONS WITH SCHRODINGER SYMMETRY 

In this section we first review, following [l^, a few of the most essential features of the 
Schrodinger symmetry. We then address the question of what can be said about the structure 
of general n-point functions just based on requiring invariance under this symmetry group, 



generalizing some of the discussion of 



15|, 



16|. 



The generators of the Schrodinger group are the number operator (or the 'mass' operator) 
M, the dilatation operator D which generates scale transformations 

D:{x,t)^{\x,\H), (2.1) 

the momentum and energy operators Pi and H which generate space and time translations 

Pi-.x^x + d, H:T^t + a, (2.2) 

the angular momentum operators Mij which generate spatial rotations 

Mij-.x^Rx, (2.3) 

the generators Ki of Galilean boosts 

K,:x^x-vt, (2.4) 



and finally the generator C of special conformal transformations 



C:(x,t)^(-^,-^) , (2.5) 



I + at I + at 

which can alternatively be expressed as C = TiHTi in terms of 

Ti:(x,t)-(f,^) • (2.6) 
For an operator O of definite conformal dimension A© we have [l^ 

[D,[K,,0]]=i{^o-l)[K,.0], 

(2.7) 

[A[C,0]]=z(Ao-2)[C,0]. 

Since Ki and C evidently lower the dimension of any operator we can define a primary 
operator to be one which satisfies the conditions 

[ir„O]=0, [C,O] = 0. (2.8) 

Starting with a primary operator one can build up a tower of descendants by repeated 
commutation with H and/or Pj, which always raises the dimension. 

Let us now consider what can be said in general about the structure of an ra-point corre- 
lation function 

/l„(l,...,n) = (Oi(fi,ti)---a(x„,t„)) (2.9) 

of primary operators. First we can use translation invariance to reduce the n coordinates 
(xj, ti) to {n—1) independent coordinates (x^j, tjj), where tij = ti—tj and Xy = (xi—Xj). Since 
the [n — 1) time variables are automatically invariant under the Ki we can build conformally 
invariant variables from them just as in more familiar relativistic CFTs. Specifically, scale 
invariance implies that only ratios such as tij/t^i may appear, while invariance under special 
conformal invariance allows only the familiar cross-ratios of the form tijt^i/tiktji, of which 
n(n — 3)/2 are independent [l6|. 

Now consider the (n — 1) vectors Xij, from which we can build a total of n{n — l)/2 
independent scalars Xij ■ xm- Then we use (n — 1) Galilean boosts similar to the procedure 



m 



16| to reduce this number to {n — l){n — 2)/2. Actually the resulting independent 



conformally invariant variables can be parametrized as 

Vij = — = - — — H ] , I < J < n, (2.10) 



and this is the form in which they will appear naturally from the AdS calculations in the 
next section. One can easily check the Vij are indeed invariant under all of the Schrodinger 
group generators. 

So in general an n-point function will always be allowed to have arbitrary functional 
dependence on a total of — 3n + 1 Shrodinger-invariant variables. The functional depen- 
dence on the remaining, non-conformally invariant variables, can be determined by solving 
the analogue of the conformal Ward identities, which are differential equations expressing 
the constraints of the symmetry on correlation functions. 



For example it is well-known 15| that the 2-point function is completely fixed up to an 
overall constant to the form 



iM ^12 



A2{l,2) = c6A,,A2tT2^'e-^ . (2.11) 
The non-relatvistic 3-point function has been shown [l^ to be determined as 

^3(1,2.3) = n<r-"''""^'e'<*4+*4'F(„,) (2.12) 

i<j 

where F is an arbitrary function. Analagously, we find that the general form of the 4-point 
function is 

^4(1,2,3,4) = I I t./ ^''e'2*i4 2*24 2 t34^^('^^^^^^^^^^2,wi3,t^23) (2.13) 

~, ^14^23 ^13^24 

where A = Aj. It is a simple exercise to check that this satisfies the relevant conformal 
Ward identities. 



III. NON-RELATIVISTIC ADS/CFT 



The metric with Schrodinger group isometry constructed 



in js, 7] is 



dt^ 2d^dt + dx^ dr^ , 



+ 



(3.1) 



where x = Xi and i = 1,2 .. .d. Let us consider a massive scalar field in this background. 
The wave equation is 



(V^ - ml)<j) = (^r''+%{-^dr) + r'{2d^dt + + d^) - m^^ 0(r,e, 



Xi,t) = 0. (3.2) 



Because ^ is a compact direction we can use 



0(r,e,Xi,t) =e^^^</)Af(r,x„t) (3.3) 
to rewrite the equation (13. 2p on the Fourier modes as 
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^""^""Mz:^) + r\2iMdt + 9^) - ) ^^(r, x„ t) = , (3.4) 



where = + m^. In order to make the equation more symmetric we can introduce a 
new coordinate r) by 

0M(r, x„ t) = j dri e-^*^''^(r, r^, x„ t) . (3.5) 

This is the same trick that was introduced in [l^, where the free massive Schrodinger 
equation was reduced to a non-massive Klein-Gordon equation via Fourier transform with 
respect to the mass. 

If ifjir, T] = ±oo, Xi, t) ^ we can integrate by parts so that (13. 4p becomes effectively 

{^^'^'^ri^dr) + r\2d,dt + df) - rS^ V^(r, r^, x,, t) = , (3.6) 
and we can further simplify the equation by introducing the x coordinates according to 

t = y^(Xo - iXn+i) , 

V = ^iXo + iXn+i) , ^^'^"^ 

■^i Xi ■ 

Then the wave equation (13. ip becomes identical to that of a scalar field in a Euclidean 
AdS background 



)-m']i>{r,Xi)=0. (3.8) 



We will use exactly same strategy in the next section to compute the bulk-to-boundary 
and bulk-to-bulk propagators in the background (13. ip . 

A. The bulk-to-boundary propagator and two-point function 

The boundary of the background (13. ip is at r = and the generator associated with 
translations along the compact ^ direction idenfied as the mass operator M = So a 
{d + 3)-dimensional bulk theory is dual to a (d + l)-dimension boundary theory j^, " 



Near the boundary, solutions of the scalar wave equation behave like 
where 



A = l + ^ + y(l + 0' + m2 + M2 (3.10) 

is related to the scaling dimension of the source O in the boundary theory. We will use the 
usual AdS/CFT recipe 0, 3| to calculate correlation functions of O. 

In order to compute the two-point function of the boundary operators we have to calculate 
the on-shell action of a massive scalar field for a solution of the classical equation of motion 

V'^M=ml(PM (3.11) 

subject to the boundary condition 

lim</.M(r,t,f) = r^^+^^-^Mt.x) , (3.12) 

r— >0 

where A is given in flS.lOp . We used the compactness of ^ (13.31) in order not to write the 
explicit ^ dependence focusing on a Fourier mode (pM- 
The relevant solution of ( 13. lip is given by 

(j)M{r,x,t) = j dxidtiK{r,x,t;xi,ti)(f)o{xi,ti) , (3.13) 

where K{r,x,t] Xi,ti) is the bulk-to-boundary propagator for metric (13. ip 

■(M\A-l -i7rA/2 / N A ,,,2,f^_^.2 

where e is the regulator. As expected this heat kernel is the solution of (13. lip which in the 
limit r — > behaves as a delta function 

^A-(d+2)^(^^ f , t; f 1, ti) 6''{x - xi)6{t - ti) . (3.15) 

Let us now show how the bulk-to-boundary propagator (I3.14p may be derived by using a 
trick similar to the one described in the previous section. Using (13. 3p and (13. 5p . we get the 
following representation for the bulk-to-boundary propagator, 

ir(r,f,t;fi,ti) = dr]e-'^^'^K{r,r],x,t;xi,ti), (3.16) 



where K{r, r], x, t; Xi,ti) is the bulk-to-boundary propagator in Euchdean AdS space (13.81) [3|, 
[l8|,Q 

K{r,r],x,t;xi,ti) = ca ( ^ , of+ — ^ [ , / ^) ' (3-17) 

with A given by f l3.10l) and ca = , d ^^^^^ • 

The integral (13.161) is not well-defined as it stands, so we introduce the usual regulator e 
by t ^ t{l — ie) which is equivalent to taking t —>■ —it as was done in ^]. This leads to 
integral representation of K{r, x, t; xi, ti) 

r \^ 

j2 + 2{t - - te)7] +{x- xi)2 ) • ^^-^^^ 

It will become useful for computing higher point correlation functions in the next section. 
The result (I3.14p follows after performing the 1] integral here. 

Returning now to the solution (13.131) . we find that when evaluating the action of this 
solution only a boundary term contributes 

S[^] = lim / d'^xdt r"^"'^ 05^0, (3.19) 

and in this manner we find 



S[4>\ = OA / d'^Xidtid'^X2dt2 M^i,h)Oiti - ta) ( f-T^ ) ' ^^^^0o(a:l, ^2), (3.20) 
where^ 

aA = ^^ • (3.21) 

vrfr(A-(f + l)) 

Thus the boundary action correctly reproduces the two-point function in non-relativistic 
CFT OTTD 

^(1,2)= dl'! ,^ . ... ^(^1-^2) (^^) e'^^^. (3.22) 

>2 



7rtr(A-(f + 1)) V'^i-'^2 



This was also found in 



0. 



^ The careful treatment of the ?' — > Umit amounts to multiplying this result by an additional factor of 



2A-(d+2) 
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B. The bulk-to-bulk propagator 

The bulk-to-bulk propagator GM{r,x,t;ri,xi,ti) for a scalar field in background (13 .ip 
satisfies 

(r''^%i:^dr.)+r\2tMdt + d^) - K + M^)^ GM{r, x,t;n, x,,t,) 

= /+=^5(r - n)6\x - xi)6{t - ti) . (3.23) 

Here we again used the compactness of C, and wrote an equation for the corresponding 
Fourier mode, in terms of which the full bulk-to-bulk propagator is 

G(r, X, t- ri, ^1, fi, ti) = e^^^^-^^^GA^r, x, t; ri, fi, U). (3.24) 

If we introduce the coordinate rj by^ 

5(M-Mi)GM(r,f,t;ri,fi,ti)= ^ dr]dvie-'^'^+'^'^^^GM{r,v,x,t;n,7]uXuh) , (3.25) 
we then find that Gjuir, r], x, t; ri,rii,Xi,ti) satisfies 

(^''^'^-(^ir^-) + r^C^^A + dl) - m^^ GM{r, V, ^, t; n, r/i, Xi, h) 

= r'^+^6{r] - 7]i)6{r - ri)6'^{x - Xi)6{t - t^) . (3.26) 

This in turn is nothing other than the equation for bulk-to-bulk propagator in AdS, 
whose solution is Q 

GM{r, V, X, t- n, r/i, xl, h) = C'a(22;)^2Fi ( f + ^ " ^' '^') ' (3-27) 



where 



ir(A)r(A-|-|) 



(47r)('^+3)/2r(2A-rf - 1) 
2rri 



r'^ + r\ + {x — xiY + 2(?7 — rii){t — ti) 



(3.29) 



^ It might appear that we have introduced one more unwanted parameter Afi, but because 
G{r,r],x,t;ri,rii,xi,ti) is invariant under translations in rj, the integral will force M = Mi. That is 
why there is a S{AI — AIi) on the left hand side. 
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By changing variables we find that bulk-to-bulk propagator may be expressed as 

GM{r,x,t;ri,xi,ti) = 



t-hj Jc ^ M 2 ' 2 2' 2' \{t-ti)u 



2 

where y'^ = r'^ + rf + {x — Xi^ and C is R -|- 2{t-ti) ^^- 
Using 

n=0 

we have 



,F,(a,6;c;.) = J^^^M^.'^ (3.31) 



J due-'^''u-^2Fi{a,h-c-{k/uf) = 

= 27rM^-ie— V (»)n(^)n (-iMA;)^- ^ 
^ (c)„n! r(A + 2n) 

27rM^-ie-^'^^/2 ^ , 1 A A M^A;^ 

= f(A) ^' 2 + 2"' Y' ■ 

Putting everything together we find that the bulk-to-bulk propagator is given by 



(3.32) 



— -Ul 

GMir,x,t;ri,Xi,h)^KAeit-ti)e ^ *-*M ) Ja-^-i ij—^J , (3-33) 



t-t 



where J is the Bessel function and the constant is 

7rC'AM<^/222^-<^/2e-^^^/2r(A - d/2) , , 

Ka = . (3.34) 

IV. HIGHER-POINT CORRELATION FUNCTIONS 
A. Three-point functions 

Now let us calculate boundary three-point functions. We consider bulk interaction ver- 
tices of the form Li — cf)^ and L2 — The corresponding three-point functions are 
respectively 

/drd'^xdtdf 
-—^^Kl{r,C,x,t;x,,h)K*{r,C,x,t;x2,t2) Ks{r,C,x,t;xs,h) , (4.1) 

A'.il, 2, 3) = J d.Kld^KlK, . (4.2) 
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Using the representation fl3.16p of the bulk-to-boundary propagator, changing variables 
and using number conservation (which comes from the ^ integral) we find 

MW) = cJ da'd^'e-^^^-'e-^^^' f ^jf^ ( ^_ -Y' 

^ ' ' ^ y ^ J r<^+^ + 2{t - ti){-i - a') + {x - xif ) 

^2 / „ \ As 

, (4.3) 



j2 + 2{t - t2)(7 - P') + X2Y) + 2(t - t3)7 + (x - x-^f 

where Ci = c/^^c/^^c/^^. 



Using the known relativistic AdS / CFT 17| result and introducing the regulator, we can 
perform the integral to arrive at 



2 2 
^'h ^13 I MSi. ^23 ^ 



where 



and 



Ai(l, 2, 3) = ai(ti2)-^^^'^/'(t23)"^^^'^/'(ti3)-^^^'^/'e^^^^^^^^/ (4.4) 



^23.1 -I 

^ r[|(A-(rf + 2))]Mf^-^M, - -y-^/^ 

^1 2A/2-%'^r[Ai - f - i]r[A2 - f - i]r[A3 - f - 1] ^ ' ^ 



A ic2 (m - + t;i2(l + «e))^-.3/2 ^a,3,i/2^a,3,2/2 ^^-^^ 

= ^(t.3)^(t,3)S(%^, %^)$,(%^, %i - 1, A„ .M,.,2) (4.7) 

^2 

where the contour is (R -|- i-^^^) and 

1 ^2 ^2 ^2 

-i2 = J(fi^ + f^-^), (4.8) 

^ ti2 t23 tl3 

A = ^ Ai, A,,- fc = A, + A,- - Afc , (4.9) 

i 

A23 1 

c^^-y-'f^y-r (4.10, 

The beta function B{u,\) and the confluent hypergeometric $i(a;, /?, 7, x, ?/) are defined 
in flA9p and flAlOp respectively. 

Similarly, for the second type of interaction we find 



,iV£i ^ M2 



^2(1, 2, 3) = {a, - aiMiM2)(ti2)-^^^'^/'(t23)-^^^'^/'(ti3)^^"'^/'e^^^^+^^)j (4.11) 
where 

02 = ai [A2A3 + hd + 2- A)(A23,i)] (4.12) 
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and / is the same integral as (14.61) . The details of the evaluation of (14. 6 p are shown in 



12 1 where the 



an appendix. We note that this is the same integral which appeared in 
nonrelativistic three-point functions in CFT where identified with response functions in 
Martin-Siggia-Rose theory, but the integral has not been explicitly evaluated there. 



B. Four-point functions 

Finally, we consider the four-point function arising from a bulk quartic interaction (j)^ 
Ai(l,2,3,4) = j Kl{r,x,t;xi,ti)K*{r,x,t;X2,t2)K*{r,x,t;xs,t3)Ki{r,x,t;X4,ti) . (4.13) 
Performing the same tricks as above we find 

Ai(l,2,3,4) = j daida2da3e-'^"''e-'^''''e-'^"''Ii, (4.14) 



where 

drd'^xdtdfi f 



li = h 



Vr2 + {x- xi)2 + 2{t - ti){fi - ai) 



X 



r 



r2 + (x-X2)2 + 2(t-t2)(/U-a2), 

^3 / „ \ A4 



^r2 + (x-X3)2 + 2(t-t3)(M-«3)J [r^ + ix-x,y + 2it-t,)fi' ^^'^^^ 



is the same integral as in the relativistic AdS/CFT case [19|, l20| and 61 = caiCa2Ca3Ca4- 

Let us consider as an example the simple case, Ai = A2 = A3 = A4 = A, for which it is 
known 



20| that 



I.- '^t-t'^'IT. r...F.(A,A;2A;l-%±g!-Asmh^.) (4.16) 



where 



are conformal cross-ratios and 



V = , C = (4.17) 

X14X23 X13X24 



4 = (4 + 2^.^-«'j)- (4.18) 



Again introducing the regulator and making the change of variables Ui = ai+ (1 + i^), 
we have 

X-4 = 2ti4(l - ie)ui , xlj = 2tij(l - ie){uij + Vij{l + ie)) for i,j <4:, (4.19) 
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where Vij is defined in f l2.10p and t] and ( change similarly. 
Finally then we arrive at 

^ r(2A - I - 1) 27r('^+2)/2 ,,Ml44,M,4,,M,44. 

'^'(^•^•^•^) = r(2A) (4i..i3.)- ^ ^ 



in terms of 



f dui [ du2 [ ciu3e-'^i"^e-'^^^"^e-'^«"^ 

Jci Jc2 Jca 

(M3(wi2 + t^i2(l+^e)))-^ / rfz2i^i(A,A;2A,l-^^l^^ -sinh^^), (4.21) 



where the contour Cj is (R+i^^e). We see the result is consistent with the general four-point 
function ^A3\f . 

We could treat the scalar exchange diagram similarly by first performing the Fourier 
transform and then doing the integral over the AdS bulk, leading to 

.,(1. 2. 3. 4) ^ / .a.VW.,'e--V--V-Ve-^'/. . ,4.22) 

where the integral I4 is well-known from the original AdS/CFT literature 21 | 
dzd'^xdtdf] dwd'^ydrdt]' 



-Ki{z, T], X, t; a', f 1, ^1)^3(2, ^, t', P', ^3, ^3) 



G{z, T], X, t] w, T]', y, t)K2{w, r]', y, r; 7', fa, i2)^4(w, 77', f, r; /x', X4, t^) . (4.23) 

After a change of variables, the integral could be brought to the form ( 12.13P as required by 
non-relativistic conformal invariance. 



Acknowledgments 

We are grateful to K. Balasubramanian, A. Jevicki, K. Jin, C. Kalousios, S. Roy and 
especially M. Spradlin for useful conversations. This work was supported in part by the US 
Department of Energy under contract DE-FG02-91ER40688 and the US National Science 
Foundation under grant PHY-0643150 CAREER and PECASE. 



APPENDIX A: THREE-POINT FUNCTION INTEGRAL 



Here we consider the integral (14.61) 



Jca (m - w fi2(l + ^e))^l2.3/2^(;A23,l/2^^Al3,2/2 
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where the contour Cj is (R + ^^e) for e > 0. 

Let us consider the case ^23 > and ti2 > 0. Combining the two denominators involving 
u with a Feynman parameter z leads to 



2 ~^ 2 ■ 

r(^)r(^) Jo 



where 



A 



dz Z^-\l - Z)^~^ I dw e-^*^2«'^-A23,i/2 ^ 

C2 



-iMxu 



du 



12,3 I "13,2 



Ci [u + z{-w + Vi2{l + ie))] 2'+ 2 



27r 



(Ml) 



^12,3 I ^13,2 



r(^ + 



■ exp [iMiz{—w + f 12(1 + ie))] exp 



2 ' 2 ' 

which we have evaluated by using the identity 



^ Ai2,3 , Ai3,2 

' 2 



/ du 



-iMu 



[u + zy T{f) 



Noting that ^^§^ + = Ai, the original integral then becomes 



27r(Mi)^i-ie 



Ai-l„-j7rAi/2 /•! 



r(^)r(^; 



^12,3 



■^13.2 



where 



B 



(M2 + zMi 



A23. 



C2 



r(%i) 



exp 



ZTT A23£ 

' 2 2 



Putting everything together then gives 

^23,1 1 



r(^)r(^)r(%i; 



exp 



A 



^12.3 



X / dz z 2 ""^(1 — 2; 



^13.2 -1 ^ Ml , ^23.1 -I 



and the final integral can be evaluated thanks to the identity |22| 

x^-^l - x)^-\l + (3x)-Pe-^''dx = B{v, A)$i(i^, p, A + z/, -(3, -p) , 

which leads to the advertised result (14 .Op in terms of the beta function 

B{x,y) = [ t''-\l-t)y-^dt 
Jo 



(A2) 



(A3) 



/ > 0, M > 0, lm{z) > . (A4) 



dz z^-\l - z)^-^e'^^''^'^^^+''^B (A5) 



(A6) 



exp[^Mi2;t;i2(l + ie)] , (A7) 



(A8) 



(A9) 
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and the confluent hypergeometric function function 

(a) 



P, 7, X,y)=f^ (^W^^rn^n _ (^10) 



m,n=0 
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